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We describe gapped single-particle and collective excitations across a superfluid to insulator quan- 
tum phase transition of particles (bosons or fermions) in a periodic potential, with an even number 
of particles per unit cell. We demonstrate that the dynamics is controlled by a quantum impurity 
problem of a localized particle interacting with the bulk critical modes. Critical exponents are deter- 
mined by a renormalization group analysis. We discuss applications to spin oscillations of ultracold 
atoms in optical lattices, and to the electronic phases in the cuprate and related compounds. 
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A number of experiments have observed coherent spin 
oscillations of trapped, ultracold spinor bosonic atoms. 
With a superfluid ground state Q, Q, IE 0, 0, Q , these 
oscillations are well described by the classical (Gross- 
Pitaevski) equations which control the time evolution of 
the multi-component atomic condensate. Recent exper- 
iments 0, have also observed oscillations in a Mott 
insulating state obtained by placing the atoms in an op- 
tical lattice; in this situation, the spin oscillations can be 
viewed as the Rabi precession [3, |9j between sharp quan- 
tum states localized within each minimum of the optical 
lattice. 

In this paper, we investigate the connection between 
these two disparate pictures of the spin oscillations, by 
describing ground state spin correlations at the quantum 
critical point separating the superfluid and insulating 
phases. Bosonic atoms with total spin F — 1 and F = 2 
display a remark ably rich variety of superfluid and insu- 
lating phases [l(J, [ill 13 > anci have a correspond- 
ing plethora of universality classes of quantum phase 
transitions separating such phases. Superfluid-insulator 
transitions are also possible for paired fermionic atoms in 
an optical lattice. We defer a more complete classification 
of the spin dynamics at such transitions to a forthcom- 
ing paper. Here we focus on a class, with an even integer 
number of particles per unit cell, which displays non- 
trivial collective behavior induced by a strong coupling 
between the spin excitations and the critical number and 
phase fluctuations of the superfluid-insulator transition; 
the class includes both the fermionic and bosonic cases. 
The problem is mapped exactly onto a quantum 'impu- 
rity' problem, which couples a single localized spin ex- 
citation (the 'impurity') to the bulk critical modes; a 
solution using the renormalization group yields new crit- 
ical exponents and scaling functions. Our theory has 
some analogies to simpler models of the Kondo and X- 
ray edge effect in metals, and opens the way towards 
observing strongly-coupled quantum impurity physics in 



ultracold atom systems. 

Our results also apply to the superfluid-insulator tran- 
sition in electronic systems, with an even integer number 
of electrons per unit cell. Such a situation can arise in 
the cuprate or related compounds, with a periodic poten- 
tial generated spontaneously by 'stripe/checkerboard' or 
charge density wave order: recent experiments in insulat- 
ing spin ladder compounds |l5| have shown that each unit 
cell contains a pair of holes. Our results predict the fre- 
quency dependence of the electron photocmission spec- 
trum across a superfluid-insulator transition in which the 
'stripe' order is present on both sides of the transition; in 
other words, for a transition between a modulated insu- 
lator and a supersolid. The predictions of the spectrum 
are for gapped single particle excitations at positions of 
the gap maxima or minima in the Brillouin zone e.g. at 
the analog of the 'antinodal' points. 

We will study superfluid-insulator transitions at which 
the energy gap to both single-particle and spin excita- 
tions remains nonzero at the transition. The superfluid 
order parameter for the transition is then necessarily a 
spin singlet. The order parameter also carries a nonzero 
particle number, or 'charge' Q. Here, and henceforth, 
'particle' refers to either a single ultracold bosonic or 
fermionic atom, or an electron (but not a Cooper pair). 
For F ^ particles, this means that the simplest case 
has an order parameter, 'F, with Q = 2 and F = 0, cor- 
responding to the annihilation operator for a spin-singlet 
pair of particles, e.g. a Cooper pair. 

For definiteness, we will develop our results in the 
context of a Bose-Hubbard model for F / bosons in 
an optical lattice, and indicate the generalization to the 
fermionic case later. For this model, we consider the 
transition from a 'spin-singlet insulator' (SSI), a Mott 
insulator with an even number of atoms per lattice site 
and no spin order, to a 'spin-singlet condensate' (SSC), in 
which singlet pairs of bosons have condensed, but there is 
no single boson condensate. We will begin with a simple 
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mean-field theory of the Bose-Hubbard model, and then 
turn to a field theory of the critical properties of this and 
the corresponding fermionic model. 

The lattice bosons are annihilated with operators di iTn 
on lattice site i and spin projection m = —F . . . F. The 
Bose-Hubbard Hamiltonian can then be written as Ti. — 
—tT + V, where T is the kinetic energy term, 



E (4. 



h.c.) 



and V is the on-site interaction: 

V = Y / [U(n l -N) 2 + J\L l \ 2 } , 



(1) 



(2) 



where hi = J2 m a ln a m is the boson number operator on 
site i, and Lj is the total spin operator on site i. We 
have made the spin-independent part of the interaction 
explicitly symmetric around N particles per site. For 
F = 1, the final term is the most general spin-dependent 
interaction, but further terms are necessary for higher 
spin. To favor spin-singlet pairing in the ground state, 
we require N to be even, and J > 0. 

In the case when t — 0, the Hamiltonian is simply a 
sum of terms acting on a single site, containing only the 
commuting operators h and |L| 2 . The ground state on 
each site is therefore a spin singlet of N bosons. 

An appropriate mean-field Hamiltonian is TL m { = V — 
— T<s, — T$ , where V is the same on-site interaction as 
in Eq. @. T$ is the standard mean- field decoupling of 
the hopping term, generalized to the case with spin, 



(3) 



where ip m is a (c-number) constant vector, which will 
be used as a variational parameter. The remaining 
terms allow for the possibility of a spin-singlet conden- 
sate through the parameters \& and <&: 



and 



* u 'i,m u 'i, — m. 



(4) 



T* = £ (-i) J 

(i,j),m 



(5) 

where the factors of (— l) F+m are Clebsch-Gordan co- 
efficients that cause the boson operators to form spin- 
singlet pairs. 

We now use the ground state of Tt m f , which we denote 
|mf), as a variational ansatz and define 



fUWv*,*,*) = (mf|W|mf) , 



(6) 



which should be minimized by varying the three param- 
eters. If this minimum occurs for vanishing values of all 
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FIG. 1: Phase diagram calculated using mean-field theory. 
The calculation has been performed for spin F — 1 and filling 
factor N = 2. The horizontal and vertical axes give the tun- 
neling strength t and the spin-dependent interaction J, both 
in units of the spin-independent part of the interaction, U. 



three parameters, then |mf) breaks no symmetries and 
the Mott insulator is favored. A nonzero value for ip m 
at the minimum corresponds to a simple 'polar conden- 
sate' (PC), breaking spin-rotation symmetry; vanishing 
tjj m but nonzero values of \& and/or $ corresponds to a 
paired spin-singlet condensate. 

Since 7i m f contains terms (within T$) that link adja- 
cent sites, it cannot be straightforwardly diagonalized, as 
in the standard mean-field theory for the spinless Bose- 
Hubbard model. To find the phase boundaries, however, 
we need only terms up to quadratic order in the varia- 
tional parameters, which can be found using perturbation 
theory. Figure^shows the phase boundaries so obtained. 

We now turn to an analysis of the correlation functions 
of gapped spin-carrying modes across the SSI/SSC tran- 
sition. We will assume we are in a regime in which the 
lowest excitation with a nonzero spin is created by an 
elementary boson operator, ip m) with Q = 1 and F^O, 
with m = —F, . . . , F. 

General symmetry arguments can be used to deduce 
the low-energy theory coupling ^ to ij) m in the vicinity of 
the superfluid-insulator transition. For the field alone, 
we have the familiar |^| 4 field theory of the superfluid- 
insulator transition, 



d d xdr 



(m 2 



r*|tf| 



(7) 



in d spatial dimensions with co-ordinate x, imaginary 
time t, and (d+ l)-dimensional derivative d 2 — + V 2 . 
The even integer number of bosons per lattice site ensures 
that there is particle-hole symmetry in Sq, in the vicinity 
of the critical point. 

The same particle-hole symmetry applies to the gap- 
ped ipm field, but it is convenient to treat separately the 
particle- like (with Q = 1) and hole- like (with Q = — 1) 
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excitations. The long-time behavior will be governed by 
those excitations just above the gap A, for which the 'non- 
relativistic' limit can be taken and the action written 

s ip = Y.j ddxdT pL - ^r y2 + A ) p-n 



+ hi id, 



2m, 



-V + A h 



> (8) 



where p m and h m are fields describing particle and hole 
excitations of ip m ~ up m + v(— l) F+m h_ m , where u and 
v are coefficients chosen to that there is no p m h m term in 
Sj,. The even integer density constraint requires that the 
gap A be the same in the particle and hole sectors, but 
the masses m Pt h are, in general, allowed to be different. 

The coupling between '5 and ip m allows the conversion 
of a particle to a hole along with the creation of a pair: 



cfxdT 



■p m h m ^) . (9) 



Our central results for the gapped excitations at the 
superfluid-insulator transition follow from the quantum 
field theory defined by Sty +S^p+S g , and we will describe 
its properties below. 

However, before we do so, we note that the same the- 
ory applies also to the superfluid-insulator transitions of 
fermions, with an even integer number of fermions per 
site: the only change is that the single particle excita- 
tions p m , h m are fermionic. For a derivation for F = 1/2 
electrons, note that the filling N = 2 per site requires the 
underlying model to have at least two bands to prevent a 
trivial filled-band ground state. Suitable orthogonal lin- 
ear combinations of the electron annihilation operators 
in these bands yield the p m and e mm > h^, respectively, (e 
is the 2x2 antisymmetric tensor). 

We are interested here in the spectral functions of gap- 
ped excitations created by exciting a finite number of p 
and h quanta: at T = 0, these will have thresholds in the 
spectral density at integer multiples of A. As in previous 
work on spin ordering transitions [lfij . we now demon- 
strate that the structure of these threshold singularities 
reduces to the solution of an associated quantum impu- 
rity problem. 

Sty is isotropic in d + 1 space- time dimensions, and 
so, at its critical point, is invariant under scaling trans- 
formations with dynamic exponent 2=1. When such a 
rescaling transformation is applied to S^, it is clear that 
the coefficients of the dispersion (2m P! h) have scaling 
dimensions — 1, and are therefore irrelevant. The critical 
behavior is therefore given by (2m Pj h) = 0, which de- 
scribes a static impurity coupling to the bulk theory of 
Sty . 

Applying the same scaling analysis to the coupling in 
S g shows that the scaling dimension of g is (3 — d)/2, so 
that it is relevant for d < 3. Below we will address the 



case of two (spatial) dimensions, and describe an expan- 
sion in e = 3 — d for the critical exponents. 

We concentrate first on single-particle excitations with 
charge Q = 1 and spin F. On the insulating side of 
the transition, where rty > 0, the Green function for the 
ip m field, G^ , will have a pole corresponding to the sta- 
ble excitation at frequency u> = A, and a continuum for 
lu > A + ijrty, at frequencies large enough for a particle 
to produce a hole along with a pair. As the transition 
is approached, rty will become smaller while A remains 
fixed, and the continuum will start closer to the pole. 
On the superfluid side, the action described in Sty must 
be rewritten in terms of amplitude and phase modes of 
the condensate. Away from the transition, the ampli- 
tude mode remains gapped, while the phase mode is the 
gapless Goldstone mode corresponding to broken phase- 
rotation symmetry in the superfluid. Coupling between 
this gapless mode and ip m will cause the continuum in G^ 
to appear for frequencies just above A. The quasiparticle 
pole remains, however, due to the factors of momentum 
appearing in the matrix element for coupling to the Gold- 
stone mode, which cause the scattering rate to decrease 
as some power of u> — A. 

Exactly at the transition, by contrast, the coupling is 
not restricted by the Goldstone theorem and the criti- 
cal \& excitations qualitatively modify the structure of 
G^. The same is true of higher-order correlation func- 
tions, such as that for p m h m ', which is the simplest spin- 
carrying but charge-neutral combination. 

These correlation functions can be found using a renor- 
malization group (RG) calculation based on the formu- 
lation as a quantum impurity problem. Our approach 
is to perform a rescaling, which leaves Sty invariant, to 
relate the correlation function evaluated at a frequency 
just above the gap, say u> = X + 5lo, to another frequency, 
uj' = A + c~ e 8u>, where I is infinitesimal. The form of the 
coupling S g allows the scale-invariance of Sty to be used 
to relate correlators at these two frequencies. 

As in the standard RG, the second stage of the calcu- 
lation involves restoring the momentum cutoff to its orig- 
inal value. To lowest order in the coupling g (or, as will 
subsequently be shown to be equivalent, in an expansion 
in e = 3 — d), the only diagram that must be calculated 
for the renormalization of G^ is the self-energy diagram 



E?(iu,) 



(10) 



where the solid lines represent ip m and the dashed line 
"J. The constant part of the diagram Sf (iw), along with 
a lower-order 'tadpole' diagram that we have omitted, 
produces a renormalization of the gap A, but this is of no 
interest here. 

It happens that there are no diagrams renormalizing 
the coupling g at this order, although such diagrams 
do a ppe ar at higher order, and can be computed as in 
Ref. 
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Since ip m has only gapped excitations, it has no ef- 
fect on the bulk scaling, which is therefore given by 
the standard results for a complex |\I>| 4 theory. Rescal- 
ing by the standard phase space factor u — u^/Sd+i 
(where S d = 2/(T(d/2)(47r) d / 2 )), there is a fixed point 
with coupling u = e/5 + 0(e 2 ) in 3 — e dimensions, 
at which the scaling dimension of the field is given by 
[*] = l-f +e 2 /100 + O(e 3 ). 

These results for the bulk, along with the self-energy 
diagram above, when extended to two loops, lead to the 
RG flow equation for the coupling to the 'impurity': 



dg e 



,-,2 ; 



-^-^-^+00(^) s ),(ii) 



where we have defined g = g{ATT^ d+1 ^ 2 /T{{d - l)/2). 
The coupling therefore approaches a fixed-point value 
with g 2 = e/2 - (tt 2 /15 - 49/100)e 2 + C(e 3 ), so that the 
perturbative expansion at this point is indeed equivalent 
to an expansion in e. 

Finally, using the wavefunction renormalization of 
we arrive at the rescaling of the Green function G'^ , which 
obeys 



G^(A + 5uj) = e- yl G^{\ + e~ e Suj) 



(12) 



where y = l — g 2 +g 4 + (D(g 6 ), with g set to the fixed point 
of Eq. (jnj; this yields y = 1 - e/2 + (6/25 - 7r 2 /15)e 2 + 
0(e 3 ). This can be iterated to give G^(A + 5uo) ~ Su'V. 
Note that, at least to this order, y < 1, so that the quasi- 
particle pole at ui = A is replaced by a power-law thresh- 
old singularity. This exponent, determining the spectral 
density of the 'photoemission' of a single hole or particle 
at a band minimum or maximum, is one of our central 
results. 

A similar calculation applies to the two-particle thresh- 
hold singularity at ui — 2A. This is associated with the 
renormalization of the T mm >p m h m i correlation function, 
where T mm i is an arbitrary matrix. It is then also neces- 
sary to calculate the following insertion diagram (at one 
loop order): 



Hi 




(13) 



which accounts for the renormalization associated with 
bringing the particle and hole operators to the same 
space-time point. The resulting calculation is very simi- 
lar and leads to the result for the two-particle propagator 



where the exponent 7/2 depends upon whether the ma- 
trix T mm > is symmetric or antisymmetric. For bosons, a 
symmetric T has yi — 1 — 2g 2 + 0(g e ), while an antisym- 
metric T has ^2 = 1 exactly; these exponents therefore 
determine the threshold singularities for excitations with 
even and odd total spin F, respectively. For fermions, 
the same results apply, but with the roles of symmet- 
ric/antisymmetric T reversed. 



This paper has presented results for a variety of thresh- 
old singularities in the spectral functions of a system un- 
dergoing a superfluid-insulator transition, with an even 
number of particles per unit cell. This has direct ap- 
plication to experiments on ultracold atoms and to the 
cuprate compounds. It would be interesting to ex- 
tend these methods to study non-equilibrium phenom- 
ena, such as those measured in Ref. |7J, by methods ex- 
plored in recent studies of non-equilibrium quantum crit- 
icality El (which have so far been limited to 1+1 
dimensions). 
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